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The low-temperature thermal properties of dielectric crystals are governed by acoustic excitations 
with large wavelengths that are well described by plane waves. This is the Debye model, which rests 
on the assumption that the medium is an elastic continuum, holds true for acoustic wavelengths 
large on the microscopic scale fixed by the interatomic spacing, and gradually breaks down on 
approaching it. Glasses are characterized as well by universal low-temperature thermal properties, 
that are however anomalous with respect to those of the corresponding crystalline phases. Related 
universal anomalies also appear in the low-frequency vibrational density of states and, despite of 
a longstanding debate, still remain poorly understood. Using molecular dynamics simulations of 
a model monatomic glass of extremely large size, we show that in glasses the structural disorder 
undermines the Debye model in a subtle way: the elastic continuum approximation for the acoustic 
excitations breaks down abruptly on the mesoscopic, medium-range-order length-scale of about ten 
interatomic spacings, where it still works well for the corresponding crystalline systems. On this 
scale, the sound velocity shows a marked reduction with respect to the macroscopic value. This 
turns out to be closely related to the universal excess over the Debye model prediction found in 
glasses at frequencies of '^l THz in the vibrational density of states or at temperatures of '^lO K in 
the specific heat. 



I. INTRODUCTION 

Glasses are structurally disordered systems. As com- 
mon experience shows, in the macroscopic limit they sup- 
port sound waves as the corresponding crystalline materi- 
als do. In fact, averaging on a length-scale large enough, 
the details of the microscopic arrangement become essen- 
tially irrelevant. This holds true for sound waves with 
long wavelengths of at least several hundreds of nanome- 
ters, corresponding to wavenumbers, q, in the 10~^ nm^^ 
range, as those probed with light scattering techniques^. 
On further increasing q, the effect of the structural disor- 
der must at one point appear. The much larger g-scale of 
few nm~^ is also well known, since it can be accessed ex- 
perimentally with inelastic x-rayS and neutron^ scatter- 
ing techniques, and numerically with molecular dynamics 
simulations^i^iii^i^iiSiiii. These studies typically probe 
the dynamic structure factor, S{q,uj), that is the space 
and time Fourier transform of the density-density cor- 
relation function. They clearly indicate the existence in 
glasses of excitations that appear in those spectra as very 
broad peaks whose position as a function of q is charac- 
terized by a sinusoidal-like dispersion curve. Thus, these 
excitations strongly recall the acoustic modes in (poly- 
)crystalline systems up to roughly one half of the pseudo 
Brillouin zone^^, i.e. down to distances corresponding 
to the interatomic spacing. For this reason, they are of- 
ten dubbed as acoustic-like. However, as their broaden- 
ing clearly indicates, they are far from being crystal-like 
modes and correspond in fact to a complex pattern of 
atomic mot ions^i^i^i^ i ^ o i ^ ^ . 



Unfortunately, experimental and numerical studies 
leave a gap between the few nanometers scale and the 
hundreds of nanometers one that is extremely difhcult to 
access. This keeps still open a number of fundamental 
questions in the physics of glasses and, more in general, 
on the nature of the vibrational excitations in disordered 
systems, i) How does the transition for the acoustic- like 
excitations look like between the small-q Debye-like be- 
haviour and the large-g regime? In other words, how does 
it happen that reasonably well defined plane waves trans- 
form, on increasing q, into a complex pattern of atomic 
motions that mirror the structural disorder? ii) Similarly 
to crystalline systems, glasses are characterized by a uni- 
versal behaviour in some fundamental low-temperature 
observables like specific heat and thermal conductivityi^. 
In particular, in the ~10 K range the specific heat shows 
an excess with respect to the T^, Debye- model prediction 
and the thermal conductivity shows a plateau. Despite of 
a longstanding debatoHisOeiios^^^^^^^^^^ 

their origin is still poorly understood. It is however 
generally accepted that they are related to the ubiqui- 
tous existence at frequencies of ~1 THz of an excess of 
modes in the vibrational density of states, g{u!), over the 
Debye-model prediction goi^) = Sw^/wl), with lod the 
Debye frequency. This excess of modes is best visible in 
the reduced density of states, g(uj)/uj'^, where it appears 
as a broad feature known as boson pea k^^'^^ . Is then 
this universal behaviour in the low-temperature thermal 
properties or in the vibrational density of states related 
to the peculiar nature of the acoustic-like excitations in 
the low- frequency range? or, alternatively, do we have 
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to imagine that glasses are characterized by additional 
low- frequency modes? in) What is the information con- 
tent that we can extract from the high-g acoustic-like 
excitations measured in inelastic scattering experiments 
or calculated numerically? In which physical properties 
arc their crystal-like features (e.g. existence of dispersion 
curves) reflected? 

The experiments that have attempted to answer these 
questions accessing the difficult 10~^-1 nm^^ q-range 
lead to contrasting interpretations. An experiment based 
on a tunnel junction technique reported'^*' linear disper- 
sion for the transverse acoustic excitations in a silica glass 
up to frequencies 50% smaller than the boson peak po- 
sition in that glass (~1 THz). This result indicated that 
the acoustic excitations are unaffected in the frequency 
range relevant for the thermal properties in the '^lO K 
range, and seemed to exclude any acoustic contribution 
to the low-temperature anomalies in the specific heat of 
glasses^. Early inelastic x-ray scattering results^ seemed 
to confirm that scenario and showed crystal-like dispers- 
ing high-frequency longitudinal acoustic-like excitations 
with a broadening increasing quadratically with the fre- 
quency. Again, a simple linear dispersion of the longitu- 
dinal acoustic-like modes was observed at frequencies cor- 
responding to the boson peak position, thus suggesting 
a smooth transition between the macroscopic and micro- 
scopic regime. Recent and more accurate inelastic x-ray 
scattering studies have however revealed that the boson 
peak marks the energy where a qualitative change takes 
place: the longitudinal acoustic-like excitations show, be- 
low the boson peak position, a marked decrease of the 
phase velocity-^ and a broadening characterized by a re- 
markable fourth-power-law frequency dependence^i^^i^. 
A similar behaviour for the broadening of the transverse 
acoustic modes had been found in a silica glass at low 
temperature and at frequencies below the boson peak po- 
sition using a tunnel junction technique^i. Conversely, a 
recent experiment using inelastic ultraviolet light scat- 
tering to measure the longitudinal acoustic modes in a 
silica glass at room temperature reported the onset of 
this peculiar w'' regime at frequencies one order of mag- 
nitude smaller than the boson peak position^: the boson 
peak would then be not directly related to this regime. 
A complex and sometimes contradictory picture seems 
therefore to come out of the experiments performed so 
far. 

Classical molecular dynamics simulations have pro- 
vided a complementary tool to study the vibrational 
properties of glasses, starting from the pioneering inves- 
tigations of Rahman and coworkers'*. The body of re- 
sults available until now supports a scenario where the 
longitudinal acoustic-like excitations seem to be largely 
decoupled from the boson peak: they show a linear dis- 
persion and a broadening that grows quadratically with 
frequency with no special feature in the frequency region 
where the boson peak appears^i^i^i^ i^^i^^ . However, these 
studies could not provide a final answer on this issue due 
to the fact that the largest wavelengths that could be 



studied so far - fixed by the simulation box size and then 
ultimately by computer power ~ were still in the range of 
few tens of interatomic spacings. As a consequence, the 
corresponding lowest frequency acoustic-like excitations 
that were accessible lied too close to the boson peak po- 
sition to allow for definite conclusions in this supposingly 
crucial frequency range. 

We report here on a study of the vibrational dynam- 
ics of the classical Lennard- Jones (LJ) monatomic glass 
model^. The reason of our choice is simple: despite of the 
fact that this system easily crystallizes below the melting 
temperature and thus requires extremely fast and experi- 
mentally out of reach quenching rates in order to prepare 
a glass starting from the melt, it is the simplest realiza- 
tion of a structural glass at our disposal. It thus provides 
the very basic ingredients of the vibrational dynamics of 
a structural glass that in other systems might be super- 
imposed to additional more complex effects. Specifically, 
we present here results for an exceptionally large simu- 
lation box containing up to A^=10^ particles and clarify 
how the acoustic modes look like in the frequency region 
where the boson peak appears. In particular, we man- 
aged to probe acoustic excitations down to frequencies 
one order of magnitude lower than the boson peak posi- 
tion. Thanks to that, the study of their g-dependence al- 
lows us to establish that the boson peak originates from 
a deformation of the dispersion curves with respect to 
the crystal case. This boils down to a direct connection 
between the boson peak and the breakdown of the De- 
bye continuum approximation for the acoustic excitations 
that takes place on a length-scale matching that of the 
medium range order of the glas o^°i^^'^^ . 



II. RESULTS AND DISCUSSION 

Details about computer simulation methods used are 
given in Appendix[Xl We have calculated both the trans- 
verse, 5t((7, ijj), and the longitudinal, 5*1,(5,0;), dynamic 
structure factors which give information on the trans- 
verse and longitudinal acoustic-like excitations, respec- 
tively. It is important to underline that while the lat- 
ter can be obtained experimentally using scattering tech- 
niques, the former can only be studied in the frequency 
range relevant here using computer simulations. In what 
follows we will use LJ units. To make connection with 
experiments, we recall that, if we model argon using the 
LJ potential, then the temperature scale is in units of 
e— 125.2 K, the length scale in units of i7=3.405 A, and 
the time scale in units of t=2.11 ps. In Fig. 1 some rep- 
resentative Sriq, spectra are reported, including those 
corresponding to the lowest q value that we could reach 
in our simulation. 

For reference, we recall that at the studied number 
density p = N/V =1.015 (V being the simulation box 
volume) the melting temperature of the LJ system is 
Tm — 1, and the glass-transition temperature Tg ~ 0.4'^-. 
For what concerns the glass that we study, the Debye 
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FIG. 1: Transverse dynamical structure factors, Sriq,!-^), for 
a LJ glass at number density p=1.015, temperature r=10~^ 
and at the indicated q values, including the smallest one ac- 
cessible using a simulation box containing ~10^ atoms. The 
Brillouin peaks shift towards higher frequencies and show a 
clear broadening on increasing q. 



it clarifies that the Debye continuum approximation for 
the acoustic excitations breaks down at frequencies com- 
parable to the boson peak position, i.e. at frequencies 
much lower than previously expected. 

Complementary information on this issue can be found 
in Fig. 2b, where the frequency-dependence of the broad- 
ening, Tt, of the transverse acoustic-like excitations is 
presented. These data clearly show two different regimes: 
a Lu'^ regime at high-frequencies, as already known from 
previous studies and associated to the structural disorder 
of the glass^iiS,, and a previously unnoticed u!^ regime 
at low frequencies. The frequency that marks the change 
of regime is very close to, though slightly lower than, the 
boson peak position (see Fig. 4b). Moreover, the low 
frequency regime appears in the same range where the 
softening of the sound velocity shows up in Fig. 2a, thus 
indicating that these two features must be connected. 
This, however, is not surprising if one recalls that the 
Brillouin position and broadening can be related to the 



frequency and wavenumber are W£)=16.2 and qD=3.92, 
respectively; the first sharp diffraction peak is at ~ 7, 
that corresponds to an average nearest neighbours dis- 
tance of ~ 2n/q,n =0.9. The spectra reported in Fig. 1 re- 
fer then to (J- values down to ~10^ times smaller than the 
border of the pseudo-Brillouin zone located at ~ 9m/2. 

The 5't(9,w) spectra are characterized by two sym- 
metric peaks (Brillouin peaks) in addition to a sharp 
elastic peak at uj —0. The Brillouin peaks can be charac- 
terized by the position of the maximum and the broad- 
ening. This information has been obtained by fitting a 
damped harmonic oscillator model to the spectral region, 
/^(g,^), around the Brillouin peaks: 



/B(q,cj) cx 



TT{q)nUq) 



{uj^~nUq))^+uj^T^q)- 



(1) 



The parameters ^ril) and Txiq) represent the charac- 
teristic frequency and broadening (FWHM) of the Bril- 
louin peaks, respectively. The parameter ^It is utilized 
to obtain the transverse sound phase velocity, CT{q) = 
nT{q)/q, and is reported in Fig. 2a as a function of fre- 
quency. 

These data show an increase with frequency (posi- 
tive dispersion) of the sound velocity for uj > 0.8 as 
already reported for the longitudinal excitations in the 
same glass^. It is interesting to observe that below that 
frequency the macroscopic lo =0 sound velocity limit is 
not directly recovered: instead a previously unnoticed re- 
gion where the phase velocity decreases with increasing 
frequency (softening) appears. This is exactly the re- 
gion where the boson peak is found in this LJ glass (see 
Fig. 4b below). The boson peak then appears not in a fre- 
quency range where the acoustic-like excitations disperse 
linearly (constant phase velocity) , at odds with what was 
previously thought^''^", but where they experience a more 
complex dispersion behaviour. This is important because 
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FIG. 2: Phase velocity and broadening of the transverse 
acoustic-like excitations in a Lennard- Jones glass. Frequency 
dependence of phase velocity (red circles, panel a) and FWHM 
(red circles, panel b) of the transverse acoustic-like excitations 
of the studied LJ glass. The dashed lines in b) emphasize 



different regimes: 



at high-frequencies and 



at low- 



frequencies. The transition between the two regimes appears 
at the frequency where the phase velocity in a) shows a min- 
imum. 
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FIG. 3: Comparison between transverse and longitudinal 
acoustic-like excitations in a Lennard- Jones glass. Frequency 
dependence of phase velocity (blue squares, panel a) and 
broadening (blue squares, panel b, left axis) of the longitu- 
dinal acoustic-like excitations of the studied LJ glass. These 
data, similarly to those in Fig. 2, have been derived by fit- 
ting a damped harmonic oscillator model to the longitudinal 
dynamic structure factor spectra obtained from molecular dy- 
namics. In panel a) the longitudinal phase velocity data are 
compared to those calculated from the corresponding trans- 
verse ones (red circles) assuming a frequency-independent 
bulk modulus 5=59. Within error-bars, we can conclude that 
the frequency-dependence of the longitudinal phase velocity 
directly refiects that of the transverse one. In panel b) the 
longitudinal excitations broadening data are compared to the 
corresponding transverse ones (red circles, right axis): the two 
sets of data can be convincingly scaled one on top of the other. 
Within error-bars, we can again conclude that the frequency- 
dependence of the longitudinal data directly reflects that of 
the transverse ones. 



real and imaginary part of a complex self-energy, respec- 
tiveljiS^. 

A similar scenario holds for the longitudinal acoustic- 
like modes. The longitudinal sound velocity data (blue 
squares) are reported in Fig. 3a, and show again a de- 
crease at low frequency followed by a positive dispersion 
for uj > 0.8. The similarity with the behaviour of the 
transverse data suggests a common origin for both. In- 
deed, the longitudinal and transverse sound velocities in 
an isotropic elastic medium are simply related by the ex- 



pression: 



, , / B(uj) 4 „ , , 
clH = W^+34M. (2) 

where p is the mass density and B is the bulk modu- 
lus. This equation simply tells that in a glass the shear 
modulus, G = pc^, the longitudinal modulus, M = pcj^, 
and the bulk modulus B are related, and only two of 
them give independent information. The low-frequency, 
macroscopic value for the bulk modulus B{uj — > 0) =59 
can be obtained from the low-frequency data for cl{lo) 
and ct{uj), and is in good agreement with a literature 
value obtained for a slightly different system^^. Assum- 
ing a frequency independent value for B, we can then esti- 
mate the longitudinal sound velocity using the transverse 
data reported in Fig. 2a. The results of this calculation 
are shown in Fig. 3a (red circles). The good correspon- 
dence between the two sets of longitudinal velocity data 
derived directly from Sl{<1,uj) or from St{(It^) via Eq. [2] 
strongly supports a scenario where: i) the bulk modulus 
is frequency independent and ii) the frequency depen- 
dence of the longitudinal sound velocity simply derives 
from that of the transverse one. This picture is further 
reinforced by the comparison shown in Fig. 3b between 
the broadening of the longitudinal acoustic-like excita- 
tions (blue squares, left axis) and that of the transverse 
ones (red circles, right axis). It is clear here that these 
two quantities, within error-bars, can be scaled one on 
top of the other. This confirms that, within the accuracy 
of our calculation, the frequency dependence that char- 
acterizes the longitudinal acoustic-like excitations comes 
into play through the shear component of the longitudi- 
nal response, while the bulk component is a mere spec- 
tator. 

The observation of the regime at low frequency 
in the acoustic attenuation for both polarizations is 
an interesting result and is predicted by several mod- 
gjg25|2i|25j37_ rpj^g simplest one can be formulated in 
terms of Rayleigh scattering of the acoustic excitations 
from some kind of structural disorder—. In order to 
make this argument quantitative one needs to identify 
the scattering units that produce it, a difficult task that 
has lead to different conclusions^^. Other approaches in- 
clude a description of the disorder in terms of a contin- 
uum model with randomly fluctuating transverse elastic 
constant a^^i^^ , or of lattice models with disorder in the 
elastic constantsii. On this basis, it has been quite puz- 
zling not to observe this strong scattering regime in pre- 
vious molecular dynamics studies^'^'^'^ '-'^^i^"'^ . It has for 
example been argued that, in contrast with lattice based 
models, this feature would be absent in realistic glasses 
due to their intrinsic high level of frustration (large inter- 
nal stress)^. The present results solve the issue, clarifying 
that a strong scattering regime is present in fact in real- 
istic structural glasses also, in agreement with a number 
of theoreticali^iiii2ii^i21 and experimentat^'^ i^'^i^'^ re- 
sults. We believe that the reason why this regime got 
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FIG. 4: Results of a normal modes analysis of a Lennard- 
Jones glass in its inherent structures, a) g-dependence of the 
lowest frequency transverse eigenvalues of the studied LJ glass 
(circles) together with the function used to empirically de- 
scribe it (full line). The initial slope of the curve corresponds 
to the q -^Q limit of the transverse phase velocity of Fig. 2a 
(dashed line). Note that here q is only an approximate quan- 
tum number, and that it is increasingly difficult to associate 
higher frequency eigenvalues to specific q values, b) Reduced 
density of states of the studied LJ glass (circles) showing the 
boson peak at a; ~ 2. The macroscopic, Debye limit (dashed 
line) is indicated as well. The reduced density of states di- 
rectly obtained from the data in Fig. 4a via Eq.|3]is reported 
(full line) up to the maximum frequency where this analysis 
is appropriate. The agreement between the two calculations 
is good - note that there is no adjustable parameter in this 
comparison. This implies that the density of states, up to the 
boson peak, is well described by the low-g inflection observed 
in Fig. 4a. 



unnoticed before has to do with the fact that, beside the 
obvious cases of too small simulation boxes, the attention 
was usually focused on the longitudinal excitations (the 
only ones experimentally accessible) where this regime is 
difficult to study in detail even with a box as large as the 
one used here, see Fig. 3b. The present results in fact 
clarify that the signature of strong scattering appears for 
both polarizations at the same frequency which implies 
that - as the longitudinal speed of sound is larger than 
the transverse one (by a factor 2.3 in the present case) 
- the strong scattering regime appears for the longitudi- 
nal polarization at lower g's than for the transverse one: 
its observation is then definitively less favorable for the 
longitudinal polarization. 

More in general, the present results clearly show 
that the macroscopic and microscopic regimes for the 
acoustic-like excitations are connected by a crossover re- 
gion where the acoustic-like dispersion curves show a 
considerable deformation with respect to a simple lin- 
ear dependence: the sound velocities decrease with in- 
creasing frequency thus directly testifying the existence 
of an abrupt breakdown of the Debye approximation in 
glasses^. This comes together with the signature of 
strong scattering for the acoustic-like excitations across 



which reasonably well defined plane waves transform into 
a complex pattern of atomic vibrations. The connection 
to the boson peak is also clear: a softening of the trans- 
verse and longitudinal sound velocities implies an excess 
in the reduced vibrational density of states above the 
Debye level. Since this softening directly appears in the 
range where the boson peak is observed, or at the cor- 
responding temperatures where the excess in the specific 
heat universally appears in glassesi^, we can conclude 
that these anomalies must have an acoustic contribution. 
In the following we make this connection quantitative. 

In order to gain further insight into these molecular dy- 
namics results, we performed a standard normal modes 
analysis of the glass in its inherent structures, i.e. we 
have diagonalized the dynamical matrix calculated in lo- 
cal minima of the potential energy landscape and de- 
rived eigenvalues and eigenstates. The lowest frequency 
eigenstates of transverse polarization are close to be plane 
waves: it is possible to associate them with a leading q 
value-'* (see Appendix [X] for details). In Fig. 4a we report 
the corresponding eigenvalues as a function of q (circles) ; 
the reduced vibrational density of states derived from 
the complete normal modes analysis is shown in Fig. 4b 
(circles) . 

Fig. 4a corresponds to a pseudo-dispersion curve: it is 
not formally a dispersion curve but only an approximate 
one since in a disordered system q is not a good quan- 
tum number, though one can still imagine that it is a 
reasonable mean to count the modes of low enough fre- 
quency. This pseudo-dispersion curve must have a slope 
for q -^Q that corresponds to the transverse velocity. In 
fact, as shown in Fig. 4a, this slope well corresponds to 
the q -^0 limit of the transverse velocity (dashed line) 
derived from Fig. 2. On increasing g, the data in Fig. 4a 
show an early departure from linearity and a clear inflec- 
tion at w ~ 2 which is the frequency where the boson peak 
appears (Fig. 4b). This inflection has been interpreted in 
terms of the failure of the classical Born approximation 
for the description of the continuum elasticity in disor- 
dered systems^'*. It appears at the same frequency where 
the molecular dynamics results for the phase velocities in 
Fig. 2 and Fig. 3 show a softening. This is not surprising 
as the acoustic-like frequencies can indeed be regarded as 
average values over a distribution of eigenvalues related 
to the broadening of the Brillouin peaks. In other words, 
the velocity decrease observed at frequencies around the 
boson peak position in the acoustic-like excitations of 
transverse polarization (Fig. 2a) and, through the shear 
component of the longitudinal modulus, in those of lon- 
gitudinal polarization (Fig. 3a), directly reflects an in- 
flection that appears in the g-dependence of the low-g 
transverse eigenvalues. 

As far as we can consider q a reasonable quantum num- 
ber, it is easy to directly estimate the vibrational density 
of states starting from the knowledge of the dispersion 
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FIG. 5: The ratio nV /Q. is reported as a function of frequency 
for both the transverse (circles) and longitudinal (squares) po- 
larization. The frequency where this ratio equals one defines 
the loffe-Regel limit. For the transverse case, this limit falls 
very close to the boson peak position, indicated by the vertical 
dashed line. Different is the case of the longitudinal polariza- 
tion, where the loffe-Regel limit appears at a frequency ~ 3 
times higher than the boson peak position. 



frequency where the decay time of the acoustic-like exci- 
tations first matches half of the corresponding oscillation 
period, and therefore marks somehow an extreme upper 
bound in frequency for the validity of a plane-waves ap- 
proach as a starting point to describe the acoustic-like 
excitations. As shown in Fig. 5, the present data confirm 
that the loffe-Regel limit for the transverse excitations is 
located close to the boson peak positio n^^i^^ . Fig. 5 also 
shows that the loffe-Regel limit is reached at different fre- 
quencies for the longitudinal and transverse excitations, 
and that the longitudinal one shows up at a frequency 
higher than the boson peak position3i. However, it is 
important to underline that this last result seems not 
to be general: in a simulation study of a silica glass the 
loffe-Regel crossover was found to appear at the same fre- 
quency for both polarizations'^. The connection between 
the loffe-Regel limit for the transverse excitations and 
the boson peak^-iS is clarified by the discussion above, 
which shows that this limit is reached in fact in the fre- 
quency range where the continuum approximation for the 
acoustic excitations breaks down. 



III. CONCLUSIONS 



curves. In a simple plane-wave approach: 
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where L and T stand for the longitudinal and transverse 
polarization, and Cd is the Debye velocity. The approxi- 
mation in Eq. 3 is justified by the fact that the longitudi- 
nal contribution to the total vibrational density of states 
is small since it scales with c|,/(2c|)=4%. It is clear that 
in order to follow this approach we shall use the pseudo- 
dispersion curve obtained from the normal modes anal- 
ysis since in that case the low-frequency eigenvalues are 
reasonably close to being plane waves. In order to com- 
pute Eq. [21 instead of directly differentiating the disper- 
sion curve, we choose to use an empirical model to fit the 
simulation data (full line through the points in Fig. 4a) 
and then to differentiate the model function. The result 
of this calculation is shown in Fig. 4b (full line) and, as 
it can be appreciated, it well describes the boson peak. 
It is worth underlining that no additional input but the 
knowledge of the dispersion curve is required to perform 
this calculation. It is then clear that the reduced density 
of states of the studied glass, up to the boson peak, orig- 
inates from a deformation of the pseudo-dispersion curve 
of eigenstates that can still be approximately described 
using plane waves. 

The molecular dynamics results presented above also 
allow us to calculate the loffe-Regel limit for the acoustic- 
like excitations, defined as the frequency where = ttF 
for both polarizations"^'^. This limit corresponds to the 



Summing up, the present simulation results shed new 
light on the well known universal anomaly observed in the 
specific heat of glasses in the T ~ 10 K temperature range 
and related to the boson peak in the vibrational den- 
sity of states at frequencies of ~1 THz. We have shown 
that in glasses the elastic continuum approximation for 
the acoustic-like excitations breaks down abruptly on the 
mesoscopic, medium-range-order length-scale of about 
ten interatomic spacings, where it still works well for 
the corresponding crystalline systems. This breakdown is 
signaled by a deformation of the pseudo-dispersion curve 
and corresponds to a marked reduction of the sound ve- 
locity on the mesoscopic scale. This turns out to be the 
closely related to the aforementioned anomalies in the 
specific heat and vibrational density of states, that can 
be finally traced back onto elastic properties specific to 
glasses. 

Finally, in order to put the present results in some per- 
spective and in connection with experiments, it is impor- 
tant to emphasize once more that they refer to a simple 
monatomic LJ glass quenched with a cooling rate out 
of reach experimentally - in other words there is no ex- 
perimental analogue of the glass studied here. Still, we 
believe that such a simple model system has the great ad- 
vantage of clearly grasping fundamental features that - 
though observed in bits and pieces in many experiments 
- are often hidden by a number of additional effects in 
real glasses. This - we believe - is the reason why the 
vibrational properties of glasses are still a debated is- 
sue after several decades of studies and discussions. For 
example, the high-frequency vibrational dynamics of the 
LJ monatomic glass has been proven to be well described 
within the harmonic approximation^'^. However, clearly 
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real glasses are anharmonic systems - in what cases will 
anharmonicity start to play a role there? Moreover, dif- 
ferently from the investigated model, real glasses are of- 
ten characterized by the presence of intramolecular or 
optic-like modes - this issue is of course (and on pur- 
pose) completely disregarded here. All of these questions 
will require studies on further and more complex models 
in order to be fully addressed. 
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APPENDIX A: NUMERICAL METHODS 

The present numerical investigation has been per- 
formed using simulation boxes of different sizes contain- 
ing up to N ~10'' atoms, interacting via a Lennard- Jones 
potential with a cutoff Tc —2.5 and periodic boundary 
conditions. This has been realized using the large scale 
massively parallel molecular dynamics computer simula- 
tion code LAMMPS^. A standard microcanonical clas- 
sical molecular dynamics simulation, carried out at the 
constant number density p ~ 1.015 and at tempera- 
ture T = 2 in the normal liquid phase is followed by 
a fast quench {dT/dt ^4x10^) down to T= 10"^. The 
quenched glass sample is relaxed for a time (dependent 
on the sample size) sufficient to have a constant total en- 
ergy. The atomic positions, r,i(i), and velocities, Vi(i), 
have then been stored for a time (again dependent on 
the sample size) sufficient to get the desired resolution 
function. The time correlation functions required to ob- 
tain the dynamic structure factor, S'l('Z, i^), and its anal- 
ogous function for the transverse excitations, ST{q,Lo) 
have been computed as: 

5„(q,c.) = ^(^)Y dt{Uq.t)-Uq,Q))e^-\ (Al) 

where a is L or T and: 

N 

hil.t) = ^[v,(0-q]qe^'»-'-'(*\ (A2) 

i=l 
N 

jT(9,0 - ^{[v,(t).q]q-v,(t)} (A3) 

i=l 



with q = q/|q|. 

Addressing the issue of the dependence of the obtained 
results on the quenching rate is a difficult task for the 
studied system due to the fact that it easily crystallizes. 
The smallest quenching rate compatible with the long 
simulation times needed to reach the desired resolution 
and with the large sizes of the simulation boxes is about 5 
decades smaller [dT/dt ~ 4x10^"^) than the one used for 
this study. We have checked that the molecular dynamics 
results reported here are independent of the quenching 
rate in this range. 

For the glass studied here a standard normal modes 
analysis has been carried out to derive the vibrational 
density of states, ^(w), from the eigenvalues of the dy- 
namical matrix calculated in the inherent structures of 
the glass. Moreover, the pseudo-dispersion curve for the 
transverse eigenvalues of the system reported in Fig. 4a 
has been constructed following Ref.— . More in detail, for 
each of the considered simulation boxes, we collected the 
four lowest frequency degenerate eigenvalues obtained 
from the diagonalization of the dynamical matrix. These 
modes correspond to the largest wavelength standing 
waves for the simulation box. Since the transverse sound 
velocity is ~2.3 times smaller than the longitudinal one, 
these four degenerate eigenvalues are all of transverse po- 
larization. For example, the first eigenvalue has a de- 
generacy of 12 and can be associated to wavevectors of 
the (±1,0,0) family; the second eigenvalue has a degen- 
eracy of 24 and can be associated to the wavevectors of 
the (±1,±1,0) family, and so on. These eigenvalues are 
size dependent: performing this analysis on simulation 
boxes of larger and larger size (up to N=256000 parti- 
cles), we then selected transverse eigenvalues with smaller 
and smaller g's. A quite broad range of frequencies and 
g's could thus be explored. This procedure can be ex- 
pected to lead to reasonable results only as far as the 
degeneracy of the eigenvalues shows up clearly enough to 
suggest a one-to-one relation to the corresponding values 
for g, and becomes less and less reliable on decreasing 
the box size or on increasing q. In the present case, we 
find that the highest q values up to which this analysis is 
still reasonable are q ^ 0.2-0.3, corresponding to a; ~ 2, 
i.e. basically up to the boson peak position, and becomes 
less and less reliable on further increasing q. 
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